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EXECUTIVE  SUMMARY 


The  velocity  measurement  technique  investigated  in  this  program  relies  on  the 
Doppler-shift  of  a  molecular  absorption  line  of  a  naturally  occuring  gaseous  species  in  a  high 
temperature  flowfield.  The  magnitude  of  the  shift  is  recorded  as  a  normalized  fluorescence 
signal,  per^nitting  measurements  to  be  made  remotely  in  the  flowfield  with  high  spatial 
resolution.  We  focus  on  the  use  of  pulsed  laser  systems  for  the  fluorescence  excitation, 
further  providing  extreme  temporal  resolution. 

•  In  the  course  of  this  program,  the  ground  work  was  laid  for  the 
development  and  demonstration  of  Planar  Laser-Induced 
Fluorescence  (PLIF)  velocimeter  which  would  measure  a  single 
component  of  velocity  in  a  plane  within  a  high-temperature  flow 
using  a  naturally  occuring  trace  species  -  NO.  This  report 
describes  the  theory  of  Doppler-shifted  LIF  and  details  the 
sensitivity  of  the  velocity  measurement  concept  to  fluctuations  in 
flow  pressure  and  temperature.  Modeling  of  candidate 
experimental  configurations  is  presented  which  predicts 
instrument  accuracy  taking  into  account  the  above  fluctuations 
and  realistic  signal  noise  levels. 

•  Considerable  effort  was  devoted  toward  the  design,  develop¬ 
ment,  and  shake-down  of  a  high-pressure,  high-temperature 
supersonic  burner.  This  facility  is  capable  of  generating  a  wide 
range  of  flow  conditions,  nominally  in  the  vicinity  of  Mach  3, 

1500  K,  and  1  atm  ■^nd  will  constitute  the  primary  test  bed  for 
Phase  II  development  efforts.  Preliminary  characterization 
experiments  were  undertaken  in  both  combusting  and  non¬ 
combusting  modes.  PLIF  images  of  NO  distributions  in  the 
exhaust  flow  were  used  to  verify  cold  flow  operation. 

•  Experiments  were  also  conducted  to  make  measurements  of  the 
detailed  NO  lineshape  in  order  to  determine  pressure  and 
thermal  broadening  parameteri  which  are  critical  for  calculation 
of  flow  velocity  from  fluorescence  signals.  A  frequency- 
doubled,  narrow-bandwidth  (single  mode)  CW  ring-dye  laser 
was  scanned  over  transitions  of  the  NO  A-X  system  in  a  static 
cell  over  a  range  of  background  gas  pressures.  The  experimen¬ 
tal  absorption  (and  fluorescence)  lineshapes  were  fit  with 
calculated  Voigt  profiles  to  establish  the  pressure-dependent 
collisional  broadening  parameter. 

A  major  result  of  this  comprehensive  Phase  I  program  was  the  development  of  a 
detailed  model  that  will  guide  in  the  design  of  the  diagnostic  instrument  that  will  be 
constructed  and  tested  in  Phase  11. 
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1.  INTRODUCTION 


The  resurgence  of  interest  in  supersonic  propulsion  systems  (SCRAMjets)  and 
controlled  hypersonic  flight  places  new  demands  on  experimental  techniques  for  high 
temperature  flow  property  measurements.  Among  the  key  parameters  in  vehicle  performance 
assessment  during  design  and  development  is  the  vector  velocity  field  around  the  body  and 
throughout  the  propulsion  system.  Current  velocity  measurement  techniques,  such  as  hot¬ 
wire  anemometry  and  laser-Doppler  anemometry  are  generally  inapplicable  in  high- 
temperature,  high-pressure,  supersonic  flow  test  facilities  or  are  subject  to  large  systematic 
errors.  The  NO  fluorescence  approach  investigated  in  this  Phase  I  SBIR  program 
circumvents  many  of  the  inherent  limitations  in  conventional  velocity  measurement 
technology  and  is  eminently  suited  to  measurements  within  SCRAMjet  propulsion  systems. 
Furthermore,  the  fluorescence  approach  will  be  capable  of  simultaneous,  multi-point 
measurements  using  laser  sheet  illumination  and  fluorescence  detection  via  sensitive,  imaging 
detector  arrays. 

The  experimental  challenges  associated  with  obtaining  velocity  measurements  in  high 
temperature,  compressible  flowfields  are  profound.  Laser-Doppler  Anemometry  (LDA) 
techniques  have  been  developed  to  the  point  where  the  three  components  of  the  mean  and 
fluctuating  velocity  vector  at  a  point  in  the  flowfield  can  be  measured.  A  single  velocity 
component  is  measured  by  interfering  two  coherent  beams  (derived  from  a  common  laser 
source)  at  a  point  in  the  flowfield  and  detecting  the  period  of  the  scattering  signature  from  a 
particle  traversing  the  resulting  fringe  pattern. 

Although  LDA  finds  common  application  in  near-isothermal,  low  Mach  number 
flows,  it  is  not  appropriate  for  all  flows.  In  reacting  flows,  the  density  variations  associated 
with  the  temperature  variations  in  the  flowfield  bias  the  velocity  statistics  toward  the  cold 
portions  of  the  flow  (where  the  density  of  scattering  centers  is  higher)  and  is  known  to  skew 
the  measured  distribution  of  rms  fluctuations.  Furthermore,  at  the  highest  temperatures 
associated  with  air-breathing  combustion  systems  (2000  to  3000  K),  the  seed  particles 
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themselves  (typically  AI2O3  or  Ti02)  breakdown  so  that  no  LDA  signal  is  generated. 

Finally,  in  supersonic  flows  with  shock  waves,  the  relatively  large  (~  1  /xm)  seed  particles 
may  not  always  track  the  velocity  of  the  gas  due  to  the  particle  inertia. 

Hot-wire  anemometry  has  been  used  in  fluid  dynamics  research  for  decades.  As  an 
intrusive  technique,  it  is  inherently  highly  perturbing  to  the  flowfield  downstream  of  the 
probe.  In  supersonic  flows,  the  perturbation  may  manifest  itself  as  shock  waves  emanating 
from  the  probe  or  the  probe  support.  Furthermore,  the  wire  (or  film)  cannot  survive  the 
high  temperatures  associated  with  combustion.  Hence,  this  technique  is  unsuitable  for  high 
temperature  reacting  flowfields.  In  addition,  the  wire  responds  to  a  combination  of  density, 
velocity,  and  temperature  so  that  independent  determination  of  any  individual  quantity  is 
ambiguous. 

Non-intrusive,  laser-based  gas  velocity  measurement  approaches  based  on  flow 
tagging  have  been  suggested  by  several  groups.  Optical  flow  tagging  approaches  are  based 
on  marking  a  line  or  a  cross  in  the  flow  with  one  laser  beam  and  then  probing  the  marked 
line  with  a  second  laser  beam  or  sheet  at  a  known  time  later.  The  displacement  of  the 
marked  gas  during  the  known  time  delay  £wes  the  average  convection  speed  during  the  time 
between  the  two  laser  pulses.  This  basic  approach  has  been  demonstrated  using  seeded 
biacetyl,  ^  nascent  molecular  or.ygen  (RELIEF  method),^  and  water  vapor. Only  the 
water  vapor  technique  has  promise  for  high  temperature  flows.  The  biacetyJ  seed  is  not 
stable  above  room  temperature  in  the  presence  of  oxygen.  The  marking  of  the  gas  in  the 
RELIEF  technique  is  accomplished  via  Raman  excitation  of  an  excited  vibrational  level  in  the 
oxgyen  ground  electronic  state  and  this  level  does  not  persist  long  enough  at  high 
temperatures  to  make  a  measurement. 

The  water  vapor  tagging  approach  uses  a  tunable  KrF  excimer  laser  to  dissociate 
H2O  into,  among  other  products,  ground  state  OH  in  a  thin  line  along  the  laser  beam.  At 
some  later  time,  a  second  laser  is  used  to  excite  this  OH  line  in  order  to  record  its 
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displacement.  No  demonstrations  to  date  have  succeeded  in  reacting  flows  where  the 
background  nascent  OH  distribution  overwhelms  the  visibility  of  the  laser-tagged  molecules. 

It  should  be  recognized  that  the  flow  tagging  measurements  are  not  a  direct,  spatially 
resolved  measurement  of  the  gas  velocity.  Indeed,  the  measurement  yields  meaningful 
results  only  in  one-dimensional,  or  nominally  one-dimensional,  flows  where  no  directional 
ambiguity  of  the  marked  gas  is  present.  In  muiti-dimensional  or  strongly  turbulent  flows,  the 
marked  line  or  cross  in  the  flowfield  will  be  rapidly  disrupted,  rendering  interpretation  of  the 
displaced  fluid  meaningless. 

The  approach  investigated  in  our  Phase  I  program  is  a  direct  measurement  of  the 
velocity  field  at  every  spatially-resolved  detection  point  with  no  restrictions  on  the 
dimensionality  or  turbulence  level  of  the  flow.  Furthermore,  the  flow  tagging  approaches 
require  multiple  wavelengths  with  concommitant  expensive  and  complex  laser  systems, 
whereas  the  technique  developed  in  this  program  is  purposely  designed  to  use  a  single, 
conventional  pulsed-dye  laser  system. 
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2.  BACKGROUND 


2.1  Laser-Induced  NO  Fluorescence 

The  approach  for  laser-induced  NO  fluorescence  pursued  in  this  program  can  best  be 
understood  by  reference  to  Figure  1,  a  partial  energy  level  diagram  for  NO.  The  figure 
shows  schematically  the  ground  X^n  electronic  state  and  the  first  excited  electronic  state, 

In  the  ground  electronic  state,  the  first  seven  vibrational  levels  are  shown.  Each 
vibrational  level  also  contains  a  manifold  of  rotational  levels  which  are  shown  schematically 
inv"=0. 


B-3345 


Figure  1.  Schematic  Energy  Level  Diagram  Showing  Important  Processes  in  Laser- 
Induced  NO  Fluorescence 
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The  vertical  arrow  at  the  figure  left  represents  the  energy  available  from  a  laser  which 
is  tuned  to  correspond  to  a  transition  from  an  isolated  rotational  level  in  the  ground 
vibrational  level  of  the  lower  electronic  state  upward  to  an  isolated  rotational  level  of  the 
ground  vibrational  level  in  the  upper  electronic  state.  These  rotational  transitions  are 
hereafter  referred  to  as  (0,0)  band  transitions,  denoting  the  vibrational  quantum  numbers  of 
the  upper  and  lov/er  state,  respectively. 

Molecules  elevated  to  the  upper  electronic  state  may  transfer  to  nearby  rotational 
levels  or  back  to  the  ground  electronic  state  (the  latter  process  is  term  quenching)  by 
collisions  with  neighboring  molecules.  The  quenching  rate  is  denoted  by  Q.  Spontaneous 
radiative  decay  competes  with  the  collisional  decay  through  all  the  allowed  transition 
pathways  from  the  excited  state,  giving  rise  to  the  fluorescence  signal.  Transitions  from  the 
upper  electronic  state  to  the  various  vibrational  levels  in  the  lower  state  are  denoted  by 
A(o  j),  where  i  represents  the  quantum  number  of  the  lower  vibrational  level. 

Transitions  on  these  off-diagonal  bands  shift  the  fluorescence  to  longer  wavelengths 
and  permit  spectral  isolation  of  the  weak  fluorescence  signal  from  the  laser  wavelength. 
Indeed,  65  percent  of  the  fluorescence  occurs  in  the  (0,1),  (0,2)  and  (0,3)  bands^.  This 
aspect  of  our  LIF  approach  is  significant  in  that  it  permits  elastic  scattering  of  the  laser 
radiation  to  be  spectrally  filtered,  while  transmitting  the  fluorescence.  Hence,  scattering 
from  particulates,  windows,  or  surfaces  in  the  test  facility  may  be  efficiently  discriminated. 

Details  of  the  rotational  level  structure  in  the  upper  and  lower  electronic  states,  along 
with  the  convention  used  throughout  this  work  in  labelling  the  rotational  branches  is  shown  in 
Figure  2.  Fach  rotational  level  of  the  upper  state  is  spin-split  into  two  components 
labelled  N'  -1- 1/2  and  N'  -  1/2,  where  N'  refers  to  the  rotational  quantum  number.  Each 
rotational  level  in  the  lower  X^Tl  state  is  split  by  a  coupling  of  the  electronic  spin  and  orbital 
angular  momentum  and  each  spin-orbit  level  is  further  split  by  lambda-doubling,  resulting  in 
four  levels  for  each  N".  The  spin  splitting  in  the  upper  state  is  less  than  3  GHz  for  the 
highest  N"  transitions  used  in  this  work  which  is  less  than  the  transition  Doppler  width^. 
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Figure  2.  Rotational  Branches  of  NO  A-X  Transitions 


Hence,  transitions  originating  from  the  same  lower  level  and  terminating  in  each  of  the  two 
spin  components  of  the  upper  state  are  unresolved  and  are  grouped  together  with 
brackets  in  Figure  2.  The  transition  labelling  employed  here  follows  that  used  by  Earls^  for 
the  intermediate  Hunds’  case  in  the  ground  electronic  state. 
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The  number  of  fluorescence  photons  collected  in  a  single  laser  pulse  can  be  expressed 


as 


Sp  j^(co)g(&))  ^  ^  Q  •  •  Tp  (1) 

CJ 

where: 

=  laser  spectral  intensity 
=  absorption  lineshape  function 
A  =  effective  radiation  rate 
Q  =  effective  quench  rate 

/3  =  Boltzmann  population  fraction  in  absorbing  state 
Vg  s  collection  volume 
Nj  =  total  number  density  of  absorbing  species 
0  =  collection  solid  angle 
7j  s  transmissivity  of  collection  optics 
Tp  s  laser  pulse  length. 

The  first  term  may  be  understood  as  an  excitation  efficiency,  expressing  the  spectral  overlap 
of  the  laser  bandwidth  and  the  molecular  absorption  bandwidth.  The  second  term  is  also  an 
efficiency:  the  efficiency  with  which  an  excited  molecule  will  radiatively  decay.  The  third 
term  is  the  total  population  of  the  absorbing  species  in  the  collection  volume  and  the  final 
term  is  the  fraction  of  fluorescence  photons  collected. 

fhe  velocimetry  approach  investigated  in  this  work  depends  only  on  the  first 
parenthetical  term:  the  relative  overlap  of  the  laser  spectral  bandwidth  and  the  molecular 
absorption  bandwidth.  As  we  explain  below,  the  velocity  measurement  is  always  made  with 
a  ratio  of  two  fluorescence  measurements  so  that  the  remaining  common  terms  cancel. 
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2.2  Doppler-Shifted  Fluorescence 


The  essence  of  LIF  velocimetry  is  based  on  the  Doppler-shifted  absorption  frequency 
of  an  ensemble  of  atoms  or  molecules.  The  concept  is  schematically  illustrated  in  Figure  3. 
A  group  of  molecules  possessing  a  net  directed  velocity,  denoted  by  the  vector  v  is 
intersected  by  a  laser  beam  traveling  with  a  direction  and  wavelength  denoted  by  the  wave 
vector  k,  both  calculated  with  respect  to  a  stationary  reference  frame  (the  laboratory  frame). 
The  projection  of  the  gas  velocity  vector  onto  the  laser  field  wave  vector  is  given  by  the 
scalar  u. 

The  moving  gas  molecules  experience  an  optical  field  which  is  Doppler-shifted 
relative  to  an  ensemble  of  stationary  molecules.  To  first  order,  the  magnitude  of  this 
Doppler-shift  is  given  by 


Figure  3.  Schematic  Illustration  of  the  Geometry  of  the  Linear  Doppler-Shift 
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(2) 


1  ±  ^ 


where  oJq  is  the  zero-velocity  absorption  frequency  and  c  is  the  velocity  of  light.  Higher 
order  terms  due  to  the  relativistic  Doppler  effect  are  on  the  order  of  (u/c)^  and  may  be 
neglected  for  practical  terrestrial  velocities. 


The  most  direct  measure  of  the  magnitude  of  the  Doppler-shift  is  made  by  scanning  a 
laser  wavelength  across  an  absorption  lineshape  and  determining  the  shift  of  the  lineshape 

Q 

with  respect  to  a  stationary  lineshape.  This  approach  was  first  suggested  by  Measures®  and 
realized  by  Zimmerman  and  Miles.^  A  drawback  to  this  approach,  however,  is  the  time 
required  to  complete  a  stable  scan  of  the  lineshape.  Most  recently,  a  modified  ring-dye  laser 
has  been  scanned  across  OH  fluorescence  lineshapes  in  a  supersonic  reacting  flow  with  an 
effective  time  resolution  of  approximately  1  ms.^^ 

Because  the  entire  absorption  lineshape  must  be  scanned,  the  fluorescence  must  be 
collected  with  a  continuously  recording  single-point  detector.  Hence,  it  is  not  possible  to 
image  the  velocity  in  a  simultaneous,  multi-point  measurement.  Finally,  the  narrow-band 
ring  lasers  used  in  these  experiments  are  extremely  sensitive  research  instruments  requiring 
careful  setup,  constant  adjustment,  and  are  extremely  sensitive  to  vibration,  temperature 
changes,  etc.  In  our  view,  ring-dye  laser  systems  are  unlikely  to  succeed  as  a  robust  test 
facility  diagnostic  system. 

Examination  of  Eq.  (1)  reveals  that  the  scanned  approach  is.  in  principle, 
unnecessary.  Since  the  magnitude  of  the  excitation  efficiency  overlap  integral  depends  on  the 
relative  frequency  between  the  laser  spectral  intensity  lineshape  and  the  molecular  absorption 
lineshape,  it  should  be  possible  to  infer  the  relative  velocity  by  a  simple  measure  of  the 
strength  of  the  fluorescence  signal,  provided  that  the  lineshapes  of  the  laser  and  absorber  are 
accurately  known.  Of  course,  the  fluorescence  signal  must  be  normalized  to  account  for 
variations  in  absorber  number  density,  fluorescence  efficiency,  etc.  This  latter  approach  is 
the  subject  of  the  present  program.  It  offers  the  advantage  of  improved  time  resolution  and 
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the  ability  to  generalize  the  measurement  concept  to  single-shot,  multi-point  measurements 
using  laser  sheet  excitation  and  intensified  detector  array  fluorescence  imaging. 


2.3  Lineshape  Broadening  Mechanisms 


Transitions  between  real  levels  in  atoms  and  molecules  do  not  occur  at  only  discrete 
frequencies,  but  over  a  narrow  range  of  frequencies  within  the  absorption  lineshape.  The 
functional  form  of  this  lineshape  depends  on  the  surrounding  thermodynamic  state  of  the  gas. 
The  first  broadening  mechanism  we  consider  is  the  thermal  broadening  of  the  lineshape  due 
to  the  random  motion  of  the  molecules  in  the  gas.  The  translational  energy  associated  with 
the  molecules  results  in  a  Maxwellian  velocity  distribution.  Each  molecule  traveling  with  a 
given  random  velocity  experiences  a  slightly  Doppler-shifted  frequency  of  the  exciting  optical 
field.  The  ensemble  effect  of  the  Maxwellian  velocity  distribution  gives  rise  to  the  so-called 
Doppler  profile,  expressed  as:^^ 


8(a)  “ 


ln2  1 
N  11  Ao 


D 


exd-4  ln2 


(w-Wo)^ 


A(o1 


(3) 


where: 

Awj)  =  Doppler-width 
and  is  given  by 


Afajj)  =  2 


2  kT  in2 


me 


1/2 


Wg  =  frequency  of  the  absorption  line-center 


(4) 


If  Doppler-broadening  were  the  only  broadening  mechanism  of  relevance,  the  lineshape 
would  depend  on  a  single  parameter  of  the  surrounding  flowfield  -  the  temperature. 


Near  atmospheric  pressure,  however,  a  second  broadening  mechanism  is  important. 
Collisions  with  neighboring  molecules  interrupt  the  coherent  interaction  of  the  molecule  with 


10 


the  optical  field  and  perturb  the  energy  levels  of  the  molecule,  introducing  the  so-called 
Lorentz  lineshape^^ 

A(j}q 

°  : - 2 - T - 2  (5) 

+  —  Aoj^ 

where  Wg  is  the  shifted  line  center  absorption  frequency  and  is  the  collision- width.  The 
collisional  shift  arises  from  the  perturbation  of  the  isolated  molecular  energy  levels  during 
the  collision.  It  is  typically  expressed  as  a  fraction  of  the  collision  width 

&)s  =  “o 

where  /3  may  be  positive  or  negative,  depending  on  the  structure  of  the  electronic  orbitals  of 
the  colliding  species. 

The  simultaneous  action  of  both  collisional  and  Doppler-broadening  results  is 
expressed  as  a  convolution  of  the  two  profiles.  The  resultant  expression  is  the  well  known 
Voigt  profile 


This  is  usually  written  in  terms  of  the  Voigt  function 

g(c.)  =  (8) 

yTTAcjj) 

where: 
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+  (?-yr 


(9) 


and 

^  =  2v^ 
y  =  2^/in2 
a  =  \/hi2 

The  '^oigt  "a-parameter"  is  thus  a  measure  of  the  relative  contribution  of  collisional  and 
Doppler  broadening  to  the  absorption  lineshape.  Since  the  probability  for  transition 
integrated  over  the  entire  lineshape  must  remain  unchanged  in  the  presence  of  broadening 
mechanisms,  the  integral  of  the  lineshape  over  all  frequencies  is  equal  to  unity. 

While  the  Doppler-width  for  a  given  transition  is  easily  calculated,  the  collisional- 
width  and  -shift  are  more  problematic.  Physically,  the  width  is  simply  proportional  to  the 
rate  of  broadening  collisions.  As  such,  it  is  dependent  on  the  composition  of  the  surrounding 
gas,  the  overall  gas  pressure,  and  temperature.  The  pressure  dependence  is  usually  separated 
by  expressing  the  collision  width  as^^ 

where  7j  is  the  collision  broadening  parameter  for  species  i,  Pj  is  the  partial  pressure  of 
species  i,  and  the  summation  is  over  all  available  collision  partners.  The  temperature 
dependent  collisional  broadening  parameter  may  vary  with  the  particular  transition  involved, 
as  well  as  the  surrounding  gas. 
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3.  HIGH-RESOLUTION  MEASUREMENTS  OF  NO  LINESHAPE  PARAMETERS 

Few  measurements  are  available  of  the  collisional  broadening  parameters  for  NO. 

Prior  to  this  study,  the  only  broadening  measurements  available  were  a  set  of  room 
temperature  measurements  made  using  spectrally-resolved  flashlamp  absorption  measure¬ 
ments.^^  The  resolution  of  this  approach  is  not  high  so  that  the  inferred  broadening  is 
questionable  and  no  temperature-dependent  data  was  available.  As  shown  in  the  preceding 
sections,  the  velocity  measurement  technique  depends  sensitively  on  the  detailed  absorption 
lineshape,  it  is  necessary  to  accurately  measure  the  temperature-dependent  collisional 
broadening  and  shift  for  various  collisional  partners  relevant  to  high-temperature 
aerodynamics  such  as  O2,  H2O,  N2,  and  NO  self-broadening. 

The  most  accurate  approach  for  detailed  lineshape  measurements  is  based  upon  laser- 
absorption  measurements  using  a  source  whose  spectral  bandwidth  is  much  less  than  the 
inherent  linewidth  to  be  measured.  The  wavelength  of  this  source  is  tuned  across  the 
lineshape  and  the  resultant  absorption  signature  accurately  reflects  the  lineshape  at  the 
measurment  temperature,  pressure,  and  collisional  environment.  The  ability  to  generate 
tunable,  ultra-narrow  frequency  laser  radiation  at  wavelengths  shorter  than  260  nm  is  not 
generally  available  and,  in  the  interest  of  efficiency,  we  elected  to  subcontract  NO  lineshape 
measurements  to  Stanford  University,  which  has  developed  a  unique  laser  source. 

A  Coherent  699  ring  dye  laser  was  run  with  Stilbene  420  dye  and  pumped  with  6W 
(all  lines  UV)  from  a  Coherent  Innova  25/5  Ar^  laser.  Using  a  custom-designed  intra-cavity 
BBO  tuning  assembly,  3  mW  of  continuous,  single-mode  laser  radiation  was  generated  in  the 
vicinty  of  225  nm.  Tuning  was  affected  using  another  custom-designed  intra-cavity  scheme. 
This  consisted  of  removing  the  standard  etalon  and  Brewster  plate  assembly  and  inserting  a 
pair  of  Brewster-angle-cut  fused-silica  rhombs  in  front  of  the  visible  output  coupler.  The 
rhombs  were  mounted  on  a  magnesium  platform  and  driven  by  a  galvanometer  assembly  at 
its  resonant  frequency  of  2  kHz.  With  this  approach,  the  single  mode  output  of  the  laser 
(spectral  bandwidth  - 1  MHz)  could  be  scanned  continuously  over  2  to  3  cm'^.  Said 
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differently,  a  laser  field  whose  bandwidth  is  less  than  1/ 1000th  that  of  typical  NO  linewidths 
could  be  scanned  over  a  wavelength  region  several  times  the  NO  linewidth  in  less  than 
300  fis. 


This  rapid-scanning  capability  permits  fully-resolved  scans  in  transient  environments 
such  as  shock  tubes,  as  described  in  Ref.  14.  Using  this  approach,  it  is  possible  to  record 
line  broadening  parameters  over  a  wide  range  of  flow  temperatures,  pressures,  and 
compositions.  In  this  Phase  I  portion  of  the  program,  the  technique  which  will  be  used  to 
measure  relevant  broadening  parameters  was  demonstrated  in  a  static  cell  of  room 
temperature  gas  mixes,  as  shown  in  Figure  4. 

The  visible  output  of  the  ring  dye  laser  was  directed  through  a  fixed-length  etalon 
(free  spectral  range  =  2  GHz)  and  onto  a  detector.  The  transmission  through  the  etalon  was 
recorded,  providing  a  continous  measure  of  the  laser  mode  quality  and  tuning  rate  during  the 
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Figure  4.  Schematic  of  High-Resolution  NO  Lineshape  Measurements 
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laser  scan.  The  UV  output  from  the  laser  was  focussed  through  a  static  cell  containing  less 
than  1  percent  NO  and  various  background  pressures  of  N2.  A  portion  of  the  UV  beam  was 
directed  onto  a  photodetector  before  entering  the  cell  and  a  second  detector  recorded  the 
transmitted  intensity. 


The  attenuation  of  the  laser  beam  is  described  by  the  Beer-Lambert  relation 


^0 


exp  (- 


(11) 


where: 

=  spectral  absorbance  at  « 

X  =  spatial  extent  of  absorption. 

The  spectral  absorbance  is  directly  proportional  to  the  number  density  of  absorbers  and, 
more  importantly  for  our  application,  to  the  absorption  lineshape  function  derived  in 
Subsection  2.3.  A  fourth  detector  was  used  to  measure  the  laser-induced  fluorescence  from  a 
portion  of  the  beam  near  the  center  of  the  cell.  A  comparison  of  the  lineshapes  inferred 
from  both  absorption  and  fluorescence  in  th'  •  scheme  provides  a  measure  of  the  quality  of 
data  with  respect  to  saturation  or  optical  thick>.ess. 

Figure  5  shows  the  output  from  the  four  detectors  as  a  function  of  time  during  a 
series  of  typical  scans.  Each  trace  was  recorded  using  a  four-channel  Nicolet  4094 
oscilloscope  sampling  at  2  MHz.  The  top  trace  is  the  incident  laser  intensity  which  varies 
during  the  course  of  a  scan  due  to  phase-matching  variations  in  the  fixed-angle  BBO  crytal. 
The  bottom  trace  is  the  difference  between  the  output  of  the  transmitted  beam  detector  and 
the  incident  beam  detector,  which  is  proportional  to  the  absorption  lineshape.  The  third  trace 
is  the  output  from  the  fluorescence  monitoring  PMT,  also  proportional  to  the  absorption 
lineshape.  The  bottom  trace  is  the  transmission  through  the  visible  etalon  which  provides 
fiducial  marks  at  fixed  frequency  spacings  as  a  function  of  time.  (Note  that,  although  the 
scan  rate  is  not  a  linear  funtion  of  time,  the  etalon  output  provides  the  correction  factor.) 
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Figure  5.  Absorption  and  Fluorescence  Data  Traces  from  High  Resolution  Static  Cell 
Measurements 


Figure  6  is  the  reduced  absorption  lineshape  function,  normalized  to  unity  at  the  line- 
center,  of  the  isolated  P2(13)  feature  at  226.7562  nm.  The  lineshape  was  recorded  at  294  K 
with  a  cell  pressure  of  48  torr  N2  and  0.4  percent  NO.  The  experimental  data  points  are 
shown  together  with  a  least  squared-error  fit  Voigt  profile.  In  this  fit,  the  only  adjustable 
parameter  is  the  coUisional  broadening  parameter  of  Eq.  (10),  which  was  determined  to  be 
0.29  cm'^  atm'^.  This  measurement  allows  extraction  of  the  a-parameter  for  NO-N2 
broadening  at  300  K,  which  was  determined  to  be  0.31  for  this  measurement.  The  lower 
trace  is  the  residual  difference  between  the  Voigt  fit  and  the  experimental  data  and  shows 
that  the  maximum  deviations  are  less  than  2  percent. 

The  slight  asymmetry  in  the  residual  plot  is  the  result  of  neglecting  higher-order 
models  of  collision  broadened  lineshapes  which  include  asymmetric  broadening  processes. 
The  simplest  of  these  higher  order  fits  is  described  by  a  Galatry  profile.  Similar  traces  at 
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Figure  6.  NO  Absorption  Lineshape  Measurement 


pressures  up  to  0.7  atm  suggests  that  the  collisional  asymmetry  remains  a  minor  contribution 
to  the  overall  lineshape  and  can  be  neglected  for  the  present  purpose. 


By  recording  similar  traces  at  various  N2  background  pressures,  it  was  possible  to 
determine  the  frequency-shift  of  the  line-center  frequency.  Hence,  the  collisional  shift  was 
measured  to  be 

—  =  -0.303  .  (12) 
The  shift/width  ratio  is  expected  to  be  independent  of  temperature. 
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In  a  more  extensive  series  of  temperature-dependent  measurements  in  a  shock-tube 
using  dilute  NO-Ar  mixtures,  DiRosa,  et  al.,^^  report  the  following  expression  for  the  a- 
parameter  temperature  variation: 

a  =  4180  PT"^-2  , 

Because  such  complete  data  for  N2  was  not  available,  we  used  this  expression  to  evaluate  the 
a-parameter  in  all  our  analytical  modelling  along  with  the  constant  shift/width  ratio  of 
Eq.  (12). 

An  example  of  an  isolated  NO  lineshape  which  would  be  recorded  by  the  much 
broader  pulsed  dye  lasers  to  be  used  in  the  velocimeter  is  shown  in  Figure  7,  a  calculation  of 
the  Qi(33)  fluorescence  lineshape  at  1  atm  and  1500  K. 
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Figure  7.  An  NO  Voigt  Lineshape  for  the  Qi(33)  Transitions  and  Its  Derivatives 
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The  laser  lineshape  was  modelled  as  a  Gaussian  with  a  full-width  at  half-maximum  of 
0.75  cm‘^.  The  Voigt  functions  were  calculated  using  a  rational  approximation  to  the 
complex  probability  function  as  derived  by  Humlicek^^  and  checked  against  the  tabulations 
of  Armstrong.  The  magnitude  of  the  collisional  shift  at  this  pressure  is  0.05  cm'^.  The 
lower  panel  shows  the  derivative  of  the  lineshape  with  respect  to  frequency,  temperature,  and 
pressure.  The  derivatives  are  each  plotted  against  separate,  arbitrary  scales  for  visualization: 
in  fact,  the  value  of  the  temperature  derivative  is  less  than  a  few  percent  of  the  pressure  and 
frequency  derivatives. 

Inspection  of  the  Voigt  derivatives  is  instructive,  particularly  from  the  point  of  view 
of  determining  the  sensitivity  of  the  convolution  lineshape  to  temperature,  pressure,  and 
frequency  changes.  Near  the  line  center,  the  magnitude  of  the  absorption  lineshape  is 
particularly  sensitive  to  pressure:  increases  in  pressure  tend  to  depress  the  lineshape  value 
due  to  overall  broadening  of  the  lineshape.  Also  at  line  center,  the  derivative  of  the 
lineshape  with  respect  to  frequency  crosses  zero.  In  this  region  of  the  lineshape,  the 
magnitude  of  the  laser  linewidth/absorption  linewidth  overlap  (proportional  to  the  excitation 
efficiency)  is  least  sensitive  to  small  changes  in  the  relative  frequency  between  the  laser  line 
center  and  the  molecular  absorption  line  center.  Thus,  the  excitation  efficiency  is  relatively 
insensitive  to  velocity  when  the  overlap  is  concentrated  in  the  vicinity  of  the  absorption  line 
center  frequency. 

Away  from  line-center,  the  sensitivity  of  the  lineshape  to  changes  in  frequency  (and, 
hence,  Doppler-shifts  in  the  relative  frequency  due  to  directed  velocity)  increases,  reaching  a 
broad  maximum  near  the  full-width,  half-maximum  (FWHM)  of  the  absorption  lineshape. 
Coincidentally,  near  this  same  region,  the  derivatives  of  the  lineshape  with  respect  to 
temperature  and  pressure  are  near  zero.  Thus,  in  this  region  of  the  absorption  lineshape,  the 
magnitude  of  the  overlap  integral  is  most  sensitive  to  small  relative  velocity  changes  and  least 
sensitive  to  changes  in  flow  temperature  and  pressure. 
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4.  ANALYTICAL  STUDY  OF  VELOCITY  ALGORITHMS 


This  section  presents  results  from  a  detailed  analytical  study  of  the  LIF  velocimetry 
approach.  The  goal  of  the  study  was  a  quantitative  understanding  of  the  uncertainty  in  the 
inferred  velocity  due  to  pressure  and  temperature  variations  in  the  flowfield  across  the 
mea-surement  plane  as  well  as  uncertaintites  introduced  by  laser  bandwidth  fluctuations  and 
signal  noise  considerations.  Temperature  and  pressure  variations  are  particularly  important 
in  compressible  flow  applications  where  duct  cross-section  changes,  shock  waves,  or 
combustion-generated  heat  release  may  be  expected  to  significantly  modify  the  these  fields 
through  the  system.  Since  the  ultimate  goal  of  the  technique  is  multiple-point  velocity 
imaging,  it  is  important  to  examine  the  sensitivity  to  these  flowfield  variations  in  the 
measurement  plane.  Practical  experimental  considerations  such  as  laser  bandwidth  stability 
and  signal  noise  may  also  be  limitations  to  accurate  velocity  measurements,  particularly  in 
PLIF  measurements  where  the  signal  shot-noise  levels  may  be  large. 

4.1  Derivation  of  Velocity  Algorithms 

The  form  of  the  algorithm  derived  to  determine  the  flow  velocity  from  LIF 
measurements  depends  primarily  on  two  criteria: 

1.  The  dimensionality  of  the  velocity  field 

2.  The  relative  magnitude  of  the  laser  spectral  bandwidth  (Awl) 
and  the  molecular  absorption  bandwidth  (Aw^). 

We  will  begin  this  section  with  the  simplest  case  of  a  one  dimensional  flowfield,  i.e.,  vj  ^  v. 
Such  a  flowfield  might  be  found  at  the  inlet  plane  of  a  vitator-driven  tunnel,  for  example. 
Although  not  strictly  required,  it  is  useful  pedagogically  to  first  assume  that  optical  access  is 
available  both  along  the  flow  direction  (velocity  sensitive)  and  perpendicular  to  the  flow 
direction  (no  velocity  sensitivity).  This  experimental  arrangement  is  illustrated  in  Figure  8. 
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Figure  8.  Schematic  Diagram  of  Optical  Access  for  One-Dimensional  Flow 
Measurements 

The  magnitude  of  the  fluorescence  signal  determined  by  the  laser  beam  crossing  the 
flow  axis  is  denoted  by  Spj  where 

SfI  ~| 

Similarly,  the  magnitude  of  the  fluorescence  from  the  laser  beam  propagating  along  the  flow 
axis  is  given  by 

Sf2  -  I 

where  cjy  is  the  Doppler-shifted  frequency  given  by  Eq.  (2). 

In  order  to  derive  simple,  closed-form  expressions  for  the  velocity  as  a  function  of  the 
fluorescence  signal,  it  is  necessary  to  make  some  assumptions  regarding  the  relative  width  of 
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the  laser  spectral  bandwidth  and  the  molecular  absorption  lineshape,  Heuristicaliy,  it  is 
convenient  to  first  assume  that  the  laser  spectral  bandwidth  is  much  smaller  than  the 
molecular  absorption  lineshape.  This  is  schematically  illustrated  in  Figure  9,  which  shows 
the  two  functions  plotted  against  frequency.  The  magnitude  of  the  overlap  integral  is 
indicated  by  the  hatched  region. 

In  this  limit,  we  may  replace  the  laser  spectral  bandwidth  with  a  delta-function  in 
frequency  and  write 

1(0))  =  I^(o)l) 

where: 

o)L  =  laser  line-center  frequency 

and  the  normalization  of  the  delta  function  is  such  that 

I  =  I  .  (17) 


Figure  9.  Molecular  Absorption  Lineshape  and  Laser  Spectral  Bandwidth  Overlap  for  a 
Narrow-Bandwidth  Laser 
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Substituting  Eq.  (16)  into  Eqs.  (14)  and  (15)  gives 


and 


%1  =  I  g(wL) 


(18) 


S(F2)  =  I  g(coL-cOv)  - 


(19) 


Expanding  the  velocity-shifted  lineshape  function,  g(wL-Wy),  in  a  Taylor  Series  about 
and  retaining  only  the  linear  term  gives 


8(<«3L-a)v)  g(a)L) 


(20) 


Using  Eq.  (20)  to  replace  the  Doppler-shifted  lineshape  function  in  Eq.  (19)  yields  the 
following  form  for  the  fluorescence  ratio: 
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Substituting  for  the  doppler-shifted  frequency,  Wy,  and  inverting  yields  the  following 
expression  for  the  flow  velocity  as  a  function  of  the  fluorescence  ratio: 
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In  this  case,  the  velocity  algorithm  is  sensitive  to  the  inverse  of  the  slope  of  the  molecular 
absorption  lineshape  normalized  by  the  magnitude  of  the  lineshape.  As  discussed  above, 
these  quantities  are  sensitive  to  the  local  temperature  and  pressure,  depending  on  the  exact 
frequency  of  the  laser  relative  to  the  molecular  lineshape  center  frequency.  This  basic 
approach  was  used  by  Hiller,  et  al.^^‘^®  in  the  early  LIE  velocimetry  literature.  Their 
approach  was  to  seed  the  flow  with  molecular  iodine  and  probe  the  fluorescence  using  a 
frequency  modulated  Ar"^  laser. 
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In  their  work,  the  sensitivity  of  the  molecular  absorption  lineshape  proved  a  barrier  to 
quantitative  measurements.  It  was  necessary  to  use  two  separate  laser  wavelengths  (derived 
from  acousto-optic  modulation  of  the  Ar"**  frequency)  in  order  to  make  a  two-point  estimate 
of  the  slope  of  the  lineshape.  Specific  laser  excitation  frequencies  were  chosen  in  order  to 
minimize  the  pressure  and  temperature  sensitivity  of  this  two-point  slope  estimate.  Because 
the  algorithm  depended  only  on  the  linear  term  of  the  linewidth  expansion,  it  was  necessary 
to  restrict  the  range  of  velocity  shifts  such  that  the  range  of  frequencies  fell  within  the  linear 
slope  of  the  molecular  absorption  lineshape. 


The  opposite  limiting  case  is  found  when  the  laser  spectral  bandwidth  is  much  larger 
than  the  molecular  absorption  lineshape,  schematically  illustrated  in  Figure  10.  In  this  limit, 
the  molecular  absorption  lineshape  is  expressed  in  terms  of  a  delta-function 


g(w)  = 


(23) 


where: 
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Figure  10.  Molecular  Absorption  Lineshape  and  Laser  Spectral  Bandwidth  for  a  Broad- 
Bandwidth  Laser 
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By  expanding  the  laser  lineshape  in  a  Taylor  Series  and  following  the  same  procedure  as 
before,  the  following  expression  is  derived: 
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Equation  (24)  is  analogous  to  Eq.  (22)  except  that  the  normalized  lineshape  in  the 
algorithm  is  now  the  laser  lineshape  rather  than  the  molecular  absorption  lineshape.  This  is 
convenient  since  the  algorithm  is  now  independent  of  flow  temperature  and  pressure 
fluctuations  and  is  only  a  function  of  the  laser  spectral  properties.  This  approach  was  used 
by  Paul  et  al.,^^  to  measure  the  velocity  distribution  in  a  cold,  supersonic  jet  seeded  with 
NO.  In  this  case,  the  low  temperature  and  pressure  of  the  flow  collapsed  the  absorption 
lineshape  and  a  broadband  pulsed  dye  laser  was  used  to  excite  the  fluorescence. 


While  an  attractive  approach  for  low  pressure,  low  temperature  flowfields,  the 
broadband  laser  limit  is  generally  not  achievable  with  commercially  available  pulsed  dye 
lasers  in  applications  near  atmospheric  pressure.  Most  pulsed  dye  lasers  operate  with 
spectral  bandwidths  (without  intra-cavity  etalons  for  bandwidth  narrowing)  between  0.3  cm*^ 
and  0.8  cm'^  when  frequency-doubled  into  the  UV.  By  comparison,  typical  NO  linewidths 
near  atmospheric  pressure  and  1500  K  are  0.3  cm'^.  Hence,  neither  limiting  condition 
described  above  can  apply. 


In  the  intermediate  case,  illustrated  in  Figure  11,  there  is  not  a  simple,  closed  form 
expression  for  the  velocity  as  a  function  of  the  ratio  of  the  fluorescence  signals.  Rather,  we 
define  an  implicit  velocity  algorithm,  as  follows: 

rf...  _  (25) 

We  now  generalize  our  basic  approach  to  the  case  of  non-directed  flows.  As  pointed 
out  by  Hiller,^®  it  is  not  necessary  that  one  of  the  fluorescence  measurements  be  directed 
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Figure  11.  Molecular  Absorption  Lineshape  and  Laser  Spectral  Bandwidth  Overlap  for 
Equivalent  Bandwidths 


across  the  flow  with  no  Doppler-shift.  Using  counter-propagating  beams,  up-shifted  and 
down-shifted  by  the  same  magnitude,  it  is  possible  to  derive  similar  velocity  expressions. 

We  denote  the  up-shifted  values  as  Sp^  and  the  down-shifted  values  as  Sp_.  For  the  case  of 
the  narrow-bandwidth  laser,  the  form  of  the  velocity  algorithm  derived  according  to  the 
procedure  described  above  is 
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Similarly,  for  the  broad-bandwidth  laser,  the  form  of  the  algorithm  is 
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In  each  case,  it  is  only  the  premultiplier  of  the  normalized  slope  which  is  different  from  the 
expressions  derived  from  the  directed  flow  cases.  Hence,  it  is  possible  to  determine  a  single 
component  of  the  velocity  field  with  optical  access  at  an  arbitrary  angle  to  the  field.  The 
component  so  derived  will  always  be  aligned  with  the  counter-propagating  beam  paths. 

In  the  intermediate  case,  again,  no  simplification  of  the  overlap  integral  is  possible 
and  the  form  of  the  velocity  algorithm  is  given  by 

,!  _  ■  I ^(a>)g(a)-a;y) 

I  I(w)8(a;+wv)d‘«»  +  | 

In  summary,  this  section  has  shown  how  the  LIF  velocimetry  approach  may  be 
sensitive  to  the  temperature  and  pressure  of  the  flowfield  due  to  the  variations  in  the 
normalized  slope  of  the  molecular  absorption  lineshape,  depending  on  the  relative  magnitude 
of  the  laser  spectral  bandwidth  and  the  molecular  absorption  lineshape.  As  the  laser  spectral 
bandwidth  increases,  the  sensitivity  of  the  overlap  to  details  of  the  fluctuation  of  the 
molecular  absorption  bandwidth  is  diminished.  The  following  section  discusses  the 
magnitude  of  the  overlap  integral  as  a  function  of  flow  velocity  for  variable  temperature, 
pressure,  and  laser  spectral  bandwidth. 

4,2  Pulsed  Laser  Excitation  Scans 

Equation  (28)  shows  that  the  implicit  velocity  algorithm  is  the  ratio  of  a  difference 
and  a  sum  of  two  points  on  the  convolution  of  the  laser  lineshape  and  the  molecular 
absorption  lineshape.  For  a  given  velocity  shift,  two  points  on  the  convolution  lineshape  are 
subtracted  from  one  another  and  this  difference  is  normalized  by  their  sum.  Hence,  a  laser 
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excitation  scan  of  the  fluorescence  spectruin  of  a  gas  at  rest  provides  the  functional  form 
from  which  points  are  extracted  in  the  velocity  algorithm  if  the  flow  is  accelerated. 

In  order  to  identify  appropriate  transitions  to  be  used  in  the  velocimeter  and  to 
understand  the  effect  of  the  spectrally  broad  pulsed  dye  laser  on  the  absorption  lineshapes,  a 
series  of  excitation  scans  were  conducted  in  a  static  cell,  according  to  the  setup  shown  in 
Figure  12.  Tunable  UV  laser  light  was  generated  using  a  Lambda-Physik  FL  3002  dye  laser 
pumped  by  an  LPX  210i  XeCl  excimer  laser.  The  output  from  Coumarin  450  dye  was  extra¬ 
cavity  frequency  doubled  in  a  single  pass  through  an  angle-tuned,  temperature-stabilized  BBO 
crystal.  Computer  control  synchrohized  the  rotation  of  the  dye  laser  grating  with  the  crystal 
angle  so  that  continuous  scans  over  ~  3  nm  were  possible  with  pulse  energies  of  2  to  4  mJ. 
The  laser  light  was  directed  through  a  cell  (unfocussed)  in  which  various  pressures  of  air  and 
NO  could  be  introduced. 


Figure  12.  Experimental  Setup  for  Pulsed  Laser  Excitation  Scans 
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The  laser  was  operated  at  repetition  rates  between  10  and  60  Hz.  Fluorescence  was 
collected  with  a  Hamamatsu  R/955  PMT  and  averaged  over  10  to  30  laser  pulsed  using  a 
Stanford  Research  Systems  Gated  Integrator  and  Boxcar  Averager.  The  laser  scan  time  was 
adjusted  so  that  the  laser  wavelength  was  stationary  during  the  1  to  3s  integration  time 
of  the  boxcar  averager.  The  averaged  fluorescence  signal  was  recorded  on  a  laboratory 
computer. 

An  example  of  the  data  obtained  in  these  experiments  is  shown  in  Figure  13,  a 
portion  of  a  room  temperature  excitadon  scan  taken  at  P  <  1  torr  in  a  sealed  cell.  The  Pj2 
bandhead  at  227  nm  is  labelled,  along  with  portions  of  the  P2,  Q2,  and  R2  branches.  The 
linewidths  of  these  fluorescence  features  are  largely  determined  by  the  ~  0.75  cm‘^  laser 
bandwidth  since  the  inherent  Doppler  width  at  300  K  is  0.098  cm‘^.  At  higher  temperatures, 
this  portion  of  the  NO  spectrum  is  poorly  resolved  and  therefore  unsuitable  for  velocimetry 
applications. 
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Figure  13.  Room  Temperature  NO  Pulsed  Laser  Excitation  Scan 
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A  more  detailed  portion  of  the  spectrum  is  shown  in  Figure  14,  showing  the  region 
between  224.68  and  224.84  nm.  The  figure  shows  room  temperature  data  in  the  dashed  line 
along  with  a  predicted  spectrum  for  a  0.75  cm"^  bandwidth  Gaussian  laser  profile.  The 
spectral  agreement  for  the  lower  N"  values  is  excellent  because  the  computer  simulated 
spectrum  uses  observed  line  positions  where  available.  At  higher  rotational  quantum 
numbers,  the  spectral  simulation  code  predicts  the  transition  wavelength  using  a  Dunham 
coefficient  expansion.^^  The  expansion  is  accurate  to  within  about  1/3  of  the  convolution 
linewidth,  as  seen  in  the  figure.  From  this  fit,  and  a  similar  fit  taken  with  1  atm  background 
air  pressure  to  check  for  saturation  effects,  the  mean  laser  bandwidth  is  determined  to  be 
0.75  cm‘^,  in  good  agreement  with  the  dye  laser  specifications.  The  experimental  trace  also 
serves  to  verify  the  accuracy  of  the  spectral  generating  code,  allowing  spectra  to  be  simulated 
at  conditions  more  relevant  to  the  proposed  instrument. 
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Figure  14.  NO  Excitation  Scan  Detail  -  300  K 
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Figure  15  shows  a  calculated  absorption  and  excitation  spectrum  for  the  same 
wavelength  range,  now  at  1500  K  and  1  atm  air.  The  broadening  was  calculated  using  the 
a-parameter  fit  of  Eq.  (13).  The  dotted  line  is  the  actual  absorption  spectrum  and  the  solid 
line  is  the  convolution  spectrum  with  the  0.75  cm'^  bandwidth  laser.  Recall  that  the  velocity 
algorithm  is  derived  from  ratios  of  points  along  the  solid  line.  Whereas  the  Qi(33)  line  w-’.s 
very  weak  in  Figure  14,  due  to  the  relatively  low  temperature  there,  it  is  a  strong,  well- 
isolated  feature  in  Figure  15.  Extensive  survey  of  high  temperature  excitation  spectral 
simulations  revealed  that  of  the  ~  3(X)  strong  transitions  at  these  temperatures,  only  a  few 
transitions  of  the  Qj  branch  in  the  vicinity  of  N"  =  33  were  sufficiently  isolated  from  one 
another  to  ensure  that  contributions  from  neighboring  lines  did  not  complicate  the  data 
reduction.  At  wavelengths  below  224.65  nm,  (1,1)  band  transitions  overlap  the  (0,0)  band 
transitions  of  interest  here.  Hence,  the  Qi(33)  line  was  used  as  the  candidate  transition  in  all 
of  the  modelling  and  signal  strength  analyses  to  follow. 


Figure  15.  NO  Excitation  Scan  Detail  -  1500  K 
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4.3  Velocity  Algorithm  Sensitivity  Analysis 

The  form  of  the  velocity  algorithm  derived  in  Subsection  4. 1  is  an  implicit 
relationship  between  the  flow  velocity  (the  desired  experimental  parameter)  and 
sum/difference  ratios  of  lineshape  overlap  integrals  (the  measured  quantities).  Because  the 
expression  is  implicit,  it  is  not  posssible  to  derive  closed-form,  algebraic  expressions  for  the 
uncertainty  in  the  inferred  velocity  as  a  function  of  the  variation  in  the  parameters  of  the 
overlap  integrals.  Nor  is  it  possible  to  invert  Eq.  (28)  so  that  a  closed  form  expression  for 
the  velocity  uncertainty  introduced  by  signal  uncertainty  is  obtained. 

Instead,  we  have  broken  the  problem  down  into  two  components.  In  Subsec¬ 
tion  4.3.1,  we  examine  the  variation  in  the  inferred  velocity  with  variations  in  flow 
temperature,  pressure,  and  laser  bandwidth  by  calculating  the  velocity  algorithm  for  a  range 
of  parameter  values.  When  plotted  together,  the  range  of  velocities  which  give  rise  to  the 
same  value  of  the  algorithm  then  corresponds  to  the  velocity  uncertainty  introduced  by  the 
given  parameter  variation. 

In  Subsection  4.3.2,  the  uncertainty  due  to  signal  noise  is  examined  explicitly.  In 
PLIF  fluorescence  measurements,  signal  noise  is  almost  always  dominated  by  photon  shot- 
noise  carried  on  the  signal  itself.  Because  the  signal  level  varies  with  the  Doppler-shift  (and 
thus  the  velocity)  the  signal/noise  ratio  will  vary  continuously  with  velocity.  A  closed  form 
for  the  uncertainty  in  the  velocity  algorithm  due  to  the  signal  uncertainty  is  derived.  The 
velocity  algorithm  is  then  plotted  against  velocity  with  a  range  of  algorithm  values 
determined  by  the  rms  standard  deviation  due  to  the  signal-noise.  As  in  the  previous  section, 
this  plot  is  examined  in  order  to  determine  the  range  of  velocities  which  give  rise  to  a  given 
value  of  the  velocity  algorithm.  When  the  spread  in  the  velocity  algorithm  due  to  shot  noise 
is  combined  with  variations  in  the  parameters  of  Subsection  4.3.1,  a  complete  map  of  the 
velocity  uncertainty  to  signal  noise  and  property  variations  is  obtained. 
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4.3.1 


Variations  due  to  Temperature.  Pressure,  and  Laser  Bandwidth 

A  computer  routine  was  constructed  which  calculates  the  velocity  algorithm  of 
Eq.  (28)  for  Voigt  lineshapes  and  Gaussian  laser  spectral  bandwidths  as  a  function  of  the  gas 
flow  velocity.  The  parameters  of  the  implicit  velocity  equation  are  the  flow  temperature  and 
pressure  (used  to  calculate  the  Voigt  profile  for  the  given  a-parameter  dependence),  the  laser 
bandwidth,  and  the  detuning  between  the  laser  line-center  and  the  absorption  line-center  at 
zero  velocity  (hereafter  referred  to  as  the  static  detuning).  The  latter  parameter  is  the  most 
easily  controlled  in  the  experimental  design,  since  the  laser  line-center  can  be  positioned 
accurately  to  with  a  small  fraction  of  the  line- width.  The  goal  of  the  first  phase  of  the 
modelling  effort  was  to  search  for  static-detunings  which  would  result  in  minimum  sensitivity 
of  the  velocity  algorithm  to  changes  in  flow  temperature  and  pressure  over  a  given  velocity 
range. 


A  conceptual  outline  of  the  calculation  is  illustrated  in  Figure  16,  which  shows  the 
relative  lineshapes  of  the  0.75  cm'^  Gaussian  laser  line  centered  at  and  the  NO  Voigt 
absorption  line  centered  at  Wq.  A  second  absorption  lineshape  is  also  shown,  illustrating  the 
effect  of  increasing  pressure  in  broadening  the  absorption  lineshape  and  shifting  it  to  lower 
frequency.  The  pressure  shift  breaks  the  frequency  domain  symmetry  of  Eq.  (28)  and  results 
in  dramatically  different  behavior  of  the  velocity  algorithm’s  sensitivity  to  pressure  depending 
upon  whether  the  laser  static  detuning  is  to  higher  frequencies  (as  shown  in  the  figure)  or  to 
lower  frequencies. 

Figure  17  shows  the  magnitude  of  the  velocity  algorithm  plotted  as  a  function  of 
relative  flow  velocity  for  the  case  of  1500  K  air.  The  curves  constitute  a  lookup  table  for  the 
velocimeter  with  the  abcissa  being  the  experimentally  measured  parameter  and  the  ordinate 
being  the  inferred  velocity.  Two  pressure  cases  are  plotted,  0.5  and  2.0  atm,  along  with  two 
static  detunings.  The  top  two  curves  show  the  optimum  static  detuning  of 
wl  -  Wq  =  0.2  cm‘^,  where  over  the  full  dynamic  range  of  3.5  km/s,  the  factor  of  four 
variation  in  pressure  (a  typical  range  for  near  atmospheric  compressible  flows)  results  in  less 
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Figure  16.  Conceptual  Outline  of  the  Effect  of  Pressure  on  the  Velocity  Doppler-shifted 
Overlap  of  the  Laser  Lineshape  and  Absorption  Lineshape 
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than  5  percent  uncertainty  in  the  inferred  velocity.  Here  and  throughout  this  report,  the 
velocity  uncertainty  is  defined  as  the  half-width  of  the  range  of  velocities  giving  rise  to  the 
same  value  of  the  velocity  algorithm.  Also,  the  velocity  calculation  is  terminated  at  3.5 
km/s.  At  this  point,  the  smallest  signal,  Sp.,  has  fallen  to  about  2  percent  of  its  peak  value. 
Hence,  this  large  velocity  dynamic  range  requires  only  a  50-1  signal  dynamic  range,  within 
the  linear  range  of  intensified  CCD  detector  systems.  For  the  symmetric,  but  opposite 
detuning,  case  of  a  static  detuning  =  -0.2  cm‘^,  i.e.,  the  laser  is  tuned  the  same  amount  to 
the  lower  frequency  side  of  the  the  absorption  lineshape,  a  much  different  situation  obtains 
and  the  measurement  would  be  essentially  meaningless  in  the  presence  of  these  pressure 
fluctuations.  The  reason  for  this  asymmetry,  of  course,  lies  in  the  action  of  the  collisional 
shifting  and  broadening  on  the  absorption  lineshape. 

With  the  "blue"  detuning  (laser  line-center  at  a  higher  frequency  than  the  absorption 
line-center),  as  shown  in  Figure  15,  the  tendency  of  the  collisional  shift  to  move  the  value  of 
the  algorithm  at  any  velocity  further  out  on  the  convolution  lineshape  is  balanced  by  the 
overall  broadening  of  the  lineshape  so  that  the  apparent  value  of  the  algorithm  as  a  fraction 
of  the  broadening  linewidth  is  changed  little.  In  the  opposite  case,  for  the  "red"  detuning, 
the  collisional  shift  tends  to  enhance  the  broadening  of  the  convolution  lineshape  by  moving 
the  value  of  the  algorithm  nearer  to  the  line-center  of  the  convolution  profile.  As  shown  in 
Figure  7,  the  center  of  the  convolution  profile  is  less  sensitive  to  velocity  (lower  values  of 
the  algorithm  for  every  velocity)  and  more  sensitive  to  pressure  (larger  changes  in  the 
algorithm  value  with  pressure  at  a  given  velocity).  It  should  be  possible  to  find  a  "red" 
static  detuning  value  which  results  in  lower  pressure  sensitivity,  but  this  value  would  lie 
much  further  out  on  the  convolution  profile  so  that  the  pressure  shifts  back  toward  line-center 
approach  the  region  of  maximum  slope  of  the  convolution  lineshape.  In  this  scheme,  all  of 
the  fluorescence  signal  would  be  much  weaker  due  to  the  reduced  overlap  in  the  wings  of  the 
convolution  lineshape  and  the  uncertainty  in  velocity  would  be  much  higher  due  to  signal 
noise  considerations,  as  well  be  discussed  in  the  following  section. 
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It  is  important  to  note  that  in  Figure  17,  and  in  ail  following  Figures  relating  to  the 
velocity  algorithm,  the  velocity  range  shown  is  a  relative  range  about  any  reference  velocity. 
It  simply  refers  to  the  velocity  of  the  gas  relative  to  the  condition  of  the  "static"  detuning. 
This  might  often  be  zero  velocity,  but  the  laser  wavelength  could  easily  be  tuned  to  corre¬ 
spond  to  the  optimum  detuning  for  gases  flowing  at  any  velocity.  The  magnitude  of  the 
velocity  algorithm  would  then  be  used  to  measure  velocity  relative  to  this  reference  condi¬ 
tion.  As  is  evident  in  Figure  17  and  in  the  Figures  to  follow,  the  velocity  uncertainty  with 
respect  to  parametric  fluctuations  tends  to  increase  (as  a  percent  of  the  measured  value)  with 
increasing  dynamic  range.  Dynamic  range  is  here  defined  as  the  range  of  velocities  plotted, 
which  would  correspond  to  the  range  of  velocities  which  could  be  imaged  in  a  single  laser 
pulse. 


The  dependence  of  the  velocity  algorithm  on  temperature  variations  is  shown  in 
Figure  18.  The  calculations  span  the  range  between  1000  and  2000  K  at  a  fixed  flow 
pressure  of  1  atm.  As  mentioned  in  previous  sections,  the  temperature  dependance  of  the 
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NO  absorption  lineshape  is  small,  with  the  increased  broadening  of  the  Doppler  component 
balanced  by  the  decrease  in  the  collisional  component.  The  uncertainty  in  the  inferred 
velocity  is  less  than  2  percent  over  the  entire  range. 

In  a  compressible  flowfield,  the  gas  pressure  and  temperature  can  not  vary  arbitrarily. 
Coupled  variations  may  be  derived  by  considering  pressure  and  temperature  ratios  across 
shocks,  combustion-generated  heat  release,  or  isentropic  flow  variations  due  to  cross- 
sectional  area  changes  or  turbulence.  We  examined  the  latter  case  explicitly.  Taking  the 
nominal  condition  to  be  1  atm,  1500  K,  and  Mach  2  (corresponding  to  a  velocity  of 
134,0{X)  cm/s),  we  permitted  the  gas  flow  to  undergo  an  isentropic  compression  to  2  atm, 
1826  K,  Mach  1.54  (velocity  =  114,000  cm/s)  and  an  isentropic  expansion  to  0.5  atm, 

640  K,  Mach  2.44  (velocity  =  121,000  cm/s).  The  velocity  algorithm  was  calculated  for 
each  of  the  three  conditions  and  the  results  are  plotted  in  Figure  19  for  the  nominal  static 
detuning  of  0.2  cm'^. 


Figure  19.  Velocity  Algorithm  for  Isentropic  Flow  Variations 
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This  figure  shows  that  the  competing  effects  of  pressure  and  temperature  balance  one 
another  so  that  the  spread  in  the  velocity  algorithm  is  less  than  in  Figure  17  for  the  same 
range  of  pressures.  For  the  given  flow  conditions,  the  velocity  varies  by  8  percent 
(expansion)  and  15  percent  (compression),  while  the  uncertainty  at  the  nominal  condition  due 
to  the  temperature  and  pressure  variations  is  only  4  percent.  Hence,  the  uncertointy 
introduced  in  the  inferred  velocity  is  substantially  less  than  the  actual  flow  velocity  changes. 
Assuming  that  free-stream  turbulent  fluctuations  could  be  modelled  as  isentropic,  this 
suggests  that  the  velocimeter  will  be  capable  of  making  fluctuating  velocity  measurements  if 
the  velocity  fluctuations  are  greater  than  about  5  percent. 

The  final  parametric  variation  we  consider  in  this  section  is  the  effect  of  fluctuating 
laser  spectral  bandwidths.  While  the  mean  bandwidth  and  profile  may  be  accurately 
established,  there  will  be  some  residual  fluctuations  in  the  laser  bandwidth  from  shot  to  shot. 
These  fluctuations  are  not  well  characterized  by  the  manufacturer  and  we  are  aware  of  no 
published  data  on  the  bandwidth  stability  of  frequency-doubled  pulsed  dye  lasers.  As  a 
means  of  assessing  the  impact  of  such  variations,  we  calculated  the  velocity  algorithm  for 
±  5  percent  variations  in  the  nominal  laser  bandwidth  of  0.75  cm’^.  These  results  are 
plotted  in  Figure  20.  The  velocity  uncertainty  at  300,000  cm/s  is  nearly  10  percent,  the 
largest  of  the  parametric  uncertainties  investigated.  This  is  due  to  the  large  contribution  to 
the  convolution  lineshape  from  the  laser  bandwidth. 

The  5  percent  bandwidth  fluctuations  used  in  this  sample  calculation  are  larger  than 
would  be  expected  for  a  typical  pulsed  dye  laser,  so  we  may  view  this  uncertainty  as  a  worst- 
case  estimate.  This  source  of  uncertainty,  however,  can  be  removed  by  measuring  the 
instantaneous  bandwidth  of  the  laser  simultaneously  with  each  PLIF  measurement  and  using 
the  instantaneous  bandwidth  in  the  data  reduction  procedure.  The  bandwidth  of  the  visible 
laser  light  is  measured  with  a  fixed-spacing  etalon  and  the  fringe  pattern  is  recorded  using  an 
unintensified  linear  diode  array.  The  UV  bandwidth  can  be  calculated  by  multiplying  the 
visible  bandwidth  by  itself.  In  the  Phase  II  portion  of  the  program,  this  procedure  for 
determining  the  bandwidth  of  the  frequency-doubled  dye  laser  will  be  verified  by  measuring 
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Figure  20.  Velocity  Algorithm  for  Fluctuating  Laser  Bandwidth 

both  the  visible  bandwidth,  as  described  above,  and  the  UV  bandwidth  using  a  Fabry-Perot 
inteferometer  and  intensified  camera  system.  We  have  performed  similar  measurements  on  a 
Nd:YAG-pumped  dye  lasers.  With  this  approach,  the  uncertainty  displayed  in  Figure  20 
may  be  reduced  or  removed  completely. 

The  results  of  this  systematic  study  of  the  effect  of  temperature,  pressure,  and  laser 
bandwidth  on  the  instantaneous  velocity  algorithm  suggest  that  the  most  important  effect  is 
the  laser  bandwidth.  As  a  consequence,  the  Phase  II  instrument  would  include  a  means  of 
shot  to  shot  measurement  of  the  laser  bandwidth.  The  following  section  considers  the  effect 
of  signal  noise. 
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4.3.2  Signal  Noise  Considerations. 


Recall  the  equation  for  the  velocity  algorithm 


Algebraically,  this  equation  may  be  manipulated  into  the  following  form: 

^(v)  = 


2Sf^ 


-  1  = 


Sf± 


(29) 


(30) 


where  Sp^  =  Sp^.  +  Sp..  For  a  PLIF  experiment,  this  form  is  preferred  since  one  camera 
would  measure  only  Sp^  while  the  second  camera  would  integrate  both  Sp^  and  Sp..  The 
advantage  of  this  is  obvious  if  one  considers  that  the  signal  noise,  aSp,  will  be  dominated  by 
signal  shot-noise  which  is  given  by 

aSp  =  (j?Sf)^^2  (31) 

where  is  the  detector  quantum  efficiency.  Hence  the  signal-to-noise  ratio  is 

^  =  .  (32) 

ffop 

Thus  integrating  both  Sp^  and  Sp,  on  the  same  detector  avoids  the  low  signal  to  noise  ratio 
inherent  in  the  Sp.  image  alone. 


From  error  propagation  theory,  the  uncertainty  in  the  velocity  algorithm  is  related  to 
the  uncertainty  in  the  measured  quantitites  according  to  the  equation 


drp 

2  2 

+ 

aSp. 

aspi 

(33) 


Performing  the  substitution  for  the  partial  derivatives  and  manipulating  the  result 
algebraically,  we  have 
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(34) 


S^±(2S 


1/2 


The  actual  detected  signal  is  the  number  of  incident  photons,  Sp,  multiplied  by  the  detector 
quantum  efficiency,  17.  Substituting  this  and  the  expression  in  Eq.  (32)  for  the  shot-noise 
limited  uncertainty  in  each  signal,  we  have  the  following  equation  for  the  velocity  algorithm 
uncertainty  due  to  the  signal  uncertainty: 


O'P  _  2 

■  (2Sp.  -  Sp±) 


(35) 


Equation  (35)  shows  that  the  fractional  uncertainty  in  the  algorithm  approaches 
infinity  as  Sp^  2Sp^,  or  as  the  velocity  approaches  zero.  The  increased  fractional 
uncertainty  here  is  due  to  the  problem  in  measuring  small  differences  between  large  numbers 
in  the  presence  of  statistical  noise. 


To  ascertain  the  impact  of  this  additional  source  of  unceratinty  on  the  inferred 
velocity,  further  simulations  of  the  velocity  algorithm  as  a  function  of  velocity  and  with  the 
parameteric  variations  discussed  above  were  undertaken.  Signal  levels  were  predicted  from 
the  formula  introduced  in  Eq.  (1). 


Calculations  of  expected  signal  levels  were  made  assuming  a  5  cm  x  6.7  cm  imaged 
region,  5  percent  NO  in  air  at  1500  K  and  1  atm,  and  typical  values  for  the  remaining 
parameters  taken  from  our  prior  experience  with  NO  PLIF  imaging  using  frequency-doubled 
excimer-pumped  dye  laser  excitation  (cf..  Ref.  23).  It  immediately  became  apparent  that  the 
signal  uncertainty  overwhelmed  the  sources  of  uncertainty  described  in  the  previous  section 
by  a  large  measure.  In  order  to  achieve  reasonable  uncertainty  levels  with  realistic 
parameters  of  Eq.  (1)  it  would  be  necessary  to  average  neighboring  pixels  on  the  detector 
array.  This  decreases  signal  uncertainty  at  the  expense  of  spatial  resolution. 
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Figure  21  is  a  calculation  of  the  velocity  algorithm  spread  for  ±  half  the  variation  of 
Eq.  (28)  when  the  detector  pixels  are  averaged  in  blocks  of  10.  For  the  intensifed  CCD 
arrays  available  at  PSI,  with  512  x  240  elements,  this  would  correspond  to  an  effective 
resolution  of  ~  100  x  100.  This  lower  resolution  is  still  useful,  representing  over  10,000 
individual  measurement  locations  in  every  laser  shot.  With  this  averaging,  the  instantaneous 
velocity  uncertainty  due  to  the  signal  noise  is  less  than  S  percent  at  2  km/s. 

A  similar  calculation  for  the  isentropic  flow  variations  of  the  previous  section  is 
shown  in  Figure  22.  Here,  the  instantaneous  velocity  uncertainty  is  on  the  order  of 
10  percent,  as  opposed  to  4  percent  in  the  earlier  calculation.  Hence,  it  may  not  be  possible 
to  image  free-stream  turbulent  fluctuations  because  of  the  inherently  low  signal  levels 
associated  with  imaging  experiments.  The  uncertainties  of  Figure  21  could  be  achieved  in 
point  measurements,  however,  where  the  signal  levels  are  much  higher  do  to  high  laser 
intensities  (concentrating  the  available  laser  energy  into  a  single  collection  volume)  and  larger 
optical  collection  efficiencies.  Lower  uncertainties  in  mean  velocity  distributions  may  be 


Figure  21.  Velocity  Algorithm  for  Shot-Noise  Limited  Signals 
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Figure  22.  Velocity  Algorithm  for  Isentropic  Flow  Variations  in  the  Presence  of  Signal 
Shot-Noise 

achieved  in  imaging  measurements  by  averaging  multiple  velocity  images  together.  For 
example,  averaging  10  velocity  images,  obtained  approximately  Is  apart  with  the  current  PSI 
imaging  systems,  would  result  in  a  noise  reduction  of  approximately  3,  which  is  the  residual 
uncertainty  associated  with  the  isentropic  variations,  which  averaging  will  not  remove. 

In  summary,  this  section  has  detailed  the  sensitivity  of  the  proposed  velocimeter  to 
unmeasured  variations  in  the  flowfield  properties,  laser  spectral  properties,  and  signal  noise. 
Laser  spectral  bandwidth  fluctuations  were  shown  to  be  the  largest  parameteric  source  of 
velocity  uncertainty  due  the  large  contribution  of  the  laser  spectral  bandwidth  to  the 
convolution  lineshape.  This  source  of  uncertainty  may  be  removed  by  monitoring  the 
instantaneous  laser  bandwidth  with  a  second  detector  arm,  as  described.  A  signal  shot-noise 
analysis  was  presented  for  typical  NO  PLIF  imaging  configurations  and  the  contributions  of 
the  shot-noise  to  the  velocity  uncertainty  derived.  By  averaging  together  groups  of 
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5x2  pixel  sub-arrays,  resulting  in  an  overall  array  dimension  of  102  x  120,  the  shot-noise 
contribution  to  the  velocity  uncertainty  at  the  reference  conditions  of  1500  K  and  1  atm  was 
reduced  to  approximately  5  percent.  Using  this  reduced  resolution  array  approach,  it  was 
shown  that  isentropic  compressions  and  expansions  covering  a  factor  of  4  in  pressure  around 
the  reference  condition  and  including  the  shot-noise  variations  resulted  in  instantaneous 
velocity  uncertainties  on  the  order  of  10  percent. 
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5.  DEVELOPMENT  OF  HIGH  PRESSURE,  SUPERSONIC  BURNER 


This  section  describes  a  significant  engineering  development  task  of  this  Phase  I 
program.  As  detailed  in  the  previous  sections,  the  LIF  velocimetry  technique  is  sensitive  to 
fiowfield  property  variations.  It  is  therefore  critical  that  the  technique  be  quantitatively 
examined  at  flow  temperatures  and  pressures  of  interest  in  the  final  application.  For  the 
purposes  of  this  program,  these  conditions  are  representative  of  supersonic  propulsion 
facilities:  approximately  1000  K  <  T  <  2000  K,  0.5  <  P  <  2  atm,  and 
3  X  10^  cm/s  <  V  <  3  X  10^  cm/s.  Furthermore,  the  fiowfield  within  the  facility  may 
include  large  pressure  and  temperature  variations  due  to  shock  waves  and  combustion¬ 
generated  heat  release.  A  velocity  imaging  technique  should  be  insensitive  to  these  property 
variations  since  simultaneous  measurements  may  not  be  available. 

Our  goal  in  this  portion  of  the  Phase  I  program  was  to  develop  a  laboratory-scale, 
stable  source  of  high  velocity  NO  at  these  conditions.  The  source  should  have  an  adjustable 
range  of  NO  number  density,  flow  temperature,  velocity,  and  pressure.  Furthermore,  since 
detailed  PLIF  measurements  are  envisioned,  the  flow  should  be  large  enough  for  the  imaging 
technique  to  resolve  spatial  variations  and  to  permit  comparison  to  scanned  pulsed-dye  laser 
measurements  at  multiple  single-point  fiowfield  locations.  Finally,  the  fiowfield  should  be 
amenable  to  analytical  determination  of  the  property  variations  so  that  the  experimental  data 
may  be  compared  to  theoretical  predictions. 

5.1  Burner  Overview 

Our  approach  for  this  problem  was  to  develop  a  high-pressure  hydrogen-air  burner. 
The  choice  of  hydrogen  combustion  is  driven  by  the  desire  to  generate  the  highest  possible 
plenum  temperatures  as  well  as  to  approximately  match  NASP-type  SCRAMJet  chemical 
flowfields.  The  high-temperature  exhaust  products  from  this  burner  can  be  accelerated 
through  a  round  orifice  into  the  atmosphere,  creating  an  underexpanded  round  jet.  This  type 
of  fiowfield  has  been  well  studied  (cf..  Refs.  24  and  25)  so  that  our  measurements  may  be 
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compared  to  other  experimental  and  analytical  studies.  For  burner  design  pressures  up  to 
100  psi,  the  exit  pressure  ratio  of  6.8  will  generate  flow  velocities  up  to  1  x  10^  cm/s.  The 
hot  exhaust  gases  will  contain  some  concentrations  of  NO,  depending  on  the  gas  temperature 
that  is  obtained  within  the  burner.  The  NO  density  may  be  independenJy  adjusted  by  added 
small  quantities  (nominally  a  few  percent  by  volume)  to  the  plenum  flow. 


A  schematic  diagram  of  the  burner  is  shown  in  Figure  23.  The  body  of  the  burner 
consists  of  a  4  in.  I.D.  copper  cylinder.  The  nozzle  was  machined  from  bronze,  providing  a 
compromise  between  high  themtal  conductivity  and  good  maching  properties.  The  exhaust 
exit  diameter  is  0.5  cm.  The  primary  fuel  and  air  injection  is  provided  through  the  bronze 
baseplate  of  the  burner.  The  main  injector  is  an  annulus  consisting  of  a  1  in.  stainless  steel 


Figure  23.  Schematic  of  High-Pressure  Hydrogen-Air  Burner 


outer  tube  flowing  hydrogen  with  a  5/8  in.  inner  tube  for  the  air.  The  primary  flow  is 
continuously  ignited  with  a  small,  1/4  in.  pilot  flow  of  hydrogen.  The  pilot  flame  is  ignited 
at  the  beginning  of  an  experiment  with  a  spark  electrode  fitted  through  the  burner  baseplate. 

Condensation  of  the  water  vapor  at  the  burner  walls  is  vented  through  a  1/4  in.  valve, 
also  located  in  the  baseplate.  This  valve  also  contains  a  relief  stem  in  the  event  of  an  over¬ 
pressure  within  the  burner.  Burner  pressure  is  monitored  with  a  pressure  tap  connected  to  a 
Bourdon-tube  pressure  gage.  The  gas  velocities  are  controlled  and  metered  with  a  rotameter 
panel.  Flow  delivery  pressures  are  also  monitored.  The  primary  hydrogen  flow  passes 
through  an  electronic  shut-off  valve  which  is  controlled  by  the  burner  pressure.  If  the 
pressure  suddenly  drops  below  the  operating  condition  (i.e.,  an  explosive  failure)  the 
hydrogen  flow  is  automatically  shut-off.  A  thermocouple  is  used  to  monitor  the  temperature 
of  the  nozzle/bumer  0-ring  seal,  which  can  be  safely  operated  up  to  about  200°C. 

A  schematic  diagram  of  the  underexpanded  jet  at  the  maximum  operating  pressure  is 
shown  in  Figure  24  (from  Ref.  16).  For  expansion  pressure  ratios  of  about  7  and  above,  the 
characteristic  barrel  shock  pattern  and  Mach  disk  is  observed.  The  flow  accelerates  from 
Mach  1  at  the  nozzle  exit  to  about  Mach  3.5  at  the  Mach  disk.  Radially  beyond  the  barrel 
shock  is  a  lower  velocity  supersonic  flow.  Axially  beyond  the  Mach  disk,  the  flow  is 
initially  subsonic,  thereafter  relaxing  to  the  ambient  pressure  field  through  a  series  of 
diamond  shock/expansion  waves.  Along  the  nozzle  centerline  the  temperature  and  pressure 
drop,  until  the  Mach  disk,  where  they  suddenly  increase.  Further  details  of  this  flowfield 
may  be  found  in  Ref.  25.  Figure  18  shows  that,  in  the  space  of  a  few  nozzle  diameters,  the 
flowfield  contains  a  variety  of  temperature,  velocity,  and  pressure  ranges.  Hence,  it  is  an 
ideal  test  flow  for  the  LIF  velocimetry  approach. 

The  original  injector  design  consisted  of  straight  tube  flow.  A  large  velocity  mis¬ 
match  between  the  outer  fuel  and  central  air  flows  was  planned  in  order  to  induce  a  high 
level  of  mixing  between  the  two  flows  and,  ideally,  efficient  and  even  burning.  Early  shake- 
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Figure  24.  Schematic  Diagram  of  Underexpanded  Jet  Flowfield  for  an  Expansion 
Ratio  ~7 

down  experiments,  however,  proved  that  this  was  not  the  case.  The  fuel  and  air  mixing  was 
very  poor,  such  that  unmixed  gases  flowed  out  of  the  nozzle  and  ignited  in  free  jet. 

The  primary  injector  was  modified  by  capping  the  annular  air  and  H2  tubes  with 
choked,  counter-swirled  "shower-head"  distribution  covers.  With  this  modification,  mixing 
within  the  burner  was  much  more  efficient,  as  evidenced  by  a  factor  of  ~2.5  increase  in 
burner  pressure  due  to  combustion.  The  burner  was  operated  up  to  its  design  pressure  of 
100  psi.  Pressure  fluctuations,  however,  persisted  at  a  level  of  approximately  ±  5  psi  on  the 
pressure  gage.  Since  the  Bourdon-tube  gage  does  not  have  good  high  frequency  response, 
we  can  assume  that  the  pressure  fluctuations  within  the  burner  itself  were  significantly  larger. 
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The  residual  fluctuations  may  be  due  to  incomplete  choking  of  the  H2  injector  or  coupling 
through  the  un-chocked  pressure  tap  or  pilot  flame.  We  plan  to  modify  the  burner  to  reduce 
these  instabilities. 

During  the  course  of  the  shakedown  testing,  the  injector  tube  overheated,  causing  the 
silver-soldered  distribution  cap  of  the  H2  tube  to  blow  off.  Unfortunately,  the  design  of  the 
injector  was  such  that  it  could  not  be  modified  without  replacing  entirely.  Within  the 
financial  and  temporal  constraints  of  the  Phase  I  portion  of  the  program,  such  a 
comprehensive  modification  was  not  feasible. 

We  continued  to  work  with  the  flow  by  supplying  high  pressure  air  through  both  the 
H2  and  air  injector  paths.  In  this  manner,  cold  (room  temperature  plenum)  air  flows  were 
generated  with  a  maximum  internal  pressure  of  50  psi.  With  this  pressure  ratio  of 
approximately  3,  the  flowfield  downstream  of  the  orifice  is  substantially  different  from  that 
shown  in  Figure  24.  A  schematic  diagram  is  shown  in  Figure  25  (taken  from  Ref.  26), 
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Figure  25.  Schematic  Diagram  of  Underexpanded  Jet  Flowfield  for  an  Expansion  Ratio  ~3 
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The  lower  expansion  ratio  does  not  produce  a  Mach  disk  or  barrel  shock.  Instead,  the  flow 
relaxes  to  the  back  ground  pressure  through  a  series  of  alternating  expansion  fans  and  shock 
waves.  The  repeated  pattern  along  the  jet  axis  is  typically  referred  to  as  the  diamond  shock 
pattern. 

5.2  PLIF  Imaging  of  NO  Distributions  in  Exhaust  Jet 

Scanned  Doppler-shifted  velocity  measurements  in  the  cold  supersonic  flowfleld  of 
Figure  24  were  not  feasible  due  to  the  thermodynamic  properties  of  this  jet.  Because  the 
plenum  temperature  was  only  room  temperature,  the  Mach  1.5  flow  would  cool  to  about 
200  K  and  expand  to  ~  10  psi.  At  these  conditions,  the  NO  absorption  linewidths  are 
collisionally  dominated  and  approximately  0.6  cm'^.  Hence,  the  convolution  linewidth 
obtained  by  scanning  the  pulsed  laser  bandwidth  over  this  feature  would  be  nearly  1  cm,'^ 
whereas  the  maximum  Doppler-shift  would  be  on  the  order  of  0.03  cm'^.  Furthermore,  the 
pressure-dependent  collisional  shift  is  on  the  order  of  0.2  cm*^,  more  than  six  times  the 
velocity  shift.  Since  we  had  no  way  of  accurately  measuring  the  pressure  in  the  jet  flow,  it 
would  not  be  possible  to  distinguish  the  small  Doppler-shift  in  the  presence  of  uncertain 
pressure-shifts. 

In  order  to  verify  the  gross  characteristics  of  the  flowfleld,  we  seeded  the  air-only 
flow  with  a  small  amount  of  NO  and  visualized  the  cold  exhaust  flow  using  PLIF  imaging. 
The  images  were  taken  using  the  excimer-pumped  tunable,  pulsed  dye  laser  and  a  PSI 
intensified  CCD  array  camera.  The  laser  beam  was  focused  into  a  thin  sheet  and  directed 
across  the  centerline  of  the  jet  flow.  At  right  angles  to  the  laser  sheet,  the  NO  fluorescence 
from  seven  laser  pulses  was  collected  and  digitized.  An  example  of  such  an  image  is  shown 
in  Figure  26.  The  NO  fluorescence  intensity  is  false-color  encoded  according  to  the  color 
lookup  table  shown.  The  laser  wavelength  was  tuned  to  the  partially  resolved  Qj 

(1  $  N"  5  3)  bandhead  at  226.257  nm.  The  region  viewed  is  4.5  x  3.5  cm  with  the 

nozzle  exit  at  the  bottom.  The  low  N’s  used  in  the  excitation  process  are  well  below  the 
peak  N  levels  for  the  jet  exit  temperature.  Hence,  the  signal  will  increase  as  the  gas  expands 
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Figure  26.  NO  PLIF  Image  in  Supersonic  Flowfield 

into  colder,  lower  density  regions  and  the  Boltzmann  population  fraction  shifts  to  lower  N’s. 
This  is  observed  in  Figure  26  as  a  gradually  increasing  signal  level  through  the  expansion 
fans  in  the  first  diameter  downstream  of  the  exit  plane.  Across  the  recompression  shocks, 
the  temperature  suddenly  increases  and  there  is  a  corresponding  instantaneous  drop  in  signal 
level.  Comparison  of  the  distribution  of  fluorescence  in  Figure  26  with  the  sche,  .atic 
diagram  of  the  mean  flow  structure  in  Figure  25  shows  that  the  expected  alternating  diamond 
pattern  of  low  and  high  density  is  present.  A  slight  tilt  of  the  jet  is  noticeable  as  well.  The 
burner  can  was  visually  straight  in  the  lab,  so  the  flow  angle  may  be  due  to  a  nozzle  non¬ 
uniformity. 

We  plan  to  substantially  modify  the  burner  in  support  of  another  program  during  the 
coming  months.  The  goal  of  the  modifications  is  to  reduce  the  pressure  instabilities  and 
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improve  the  uniformity  of  the  exhaust  flow.  Among  the  modifications  we  will  be 
investigating  are  increasing  the  number  of  fuel/air  mixing  ports  or  platinum  catalyst-stabilized 
combustion.  We  anticipate  that  the  burner  shakedown  and  evaluation  studies  will  be 
completed  prior  to  the  start  of  the  Phase  II  portion  of  this  program. 
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6.  CONCLUSIONS  AND  SUMMARY 


This  Phase  I  SBIR  program  successfully  demonstrated  the  feasibility  of  several  key 
elements  of  a  LIF  velocimeter  based  on  NO  fluorescence.  An  analytical  treatment  of 
Doppler-shifted  fluorescence  was  dweloped  which  expresses  the  various  approaches  for 
velocity  measurements  in  a  unified  context.  Using  this  derivation,  a  detailed  model  of  the 
Doppler-shifted  fluorescence  approach  necessary  for  single-shot  velocity  imaging  was 
constructed.  An  extensive  modelling  study  showed  that  a  single  broadband,  pulsed  dye  laser 
can  obtain  time-resolved  velocity  measurements  without  an  in-situ,  simultaneous  measure¬ 
ment  of  temperature  or  pressure.  Pressure  variations  from  0.5  to  2  atm  were  shown  to  have 
a  less  than  10  percent  effect  on  the  velocity  inferred  from  the  fluorescence  signal.  Tempera¬ 
ture  variations  from  1000  to  2000  K  had  a  negligible  impact  on  the  inferred  velocity. 

A  full  simulation  of  the  proposed  imaging  velocimeter  was  carried  out  for  typical 
isentropic  flow  property  variations  about  a  nominal  Mach  2,  1500  K,  atmospheric  pressure 
air  flow.  The  simulation  calculated  the  expected  signal  levels  based  on  realistic  imaging 
configurations,  available  laser  power,  and  known  spectroscopic  data  so  that  signal  shot-noise 
effects  could  be  included  in  the  accuracy  analysis.  For  flows  where  the  pressure  and 
temperature  field  were  known,  the  single-shot  velocity  uncertainty  over  a  3.5  km/s  dynamic 
range  was  less  than  5  percent.  For  isentropic  flow  variations  over  a  factor  of  four  pressure 
variation,  the  single-shot  uncertainty  was  on  the  order  of  10  percent.  These  accuracies  were 
for  an  instantaneous  imaging  resolution  of  ~  100  x  100  pixels  over  a  5  cm  field  of  view. 

Our  single  laser  approach  greatly  simplifies  the  hardware  necessary  to  carry  out  a 
velocity  measurement.  The  simplified  setup  increases  the  number  of  available  users  of  the 
velocimetry  technique  since  single  laser  systems  are  more  commonly  found  in  test  facility 
installations.  Of  course,  additional  lasers  may  be  integrated  with  the  single  laser  approach 
described  here  in  order  to  simultaneously  measure  temperature  and  pressure.  In  many  flows, 
however,  the  additional  complexity  and  optical  access  required  for  these  approaches  may  not 
be  feasible. 
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Several  key  experimental  results  were  also  obtained.  A  novel  frequency-doubled  ring- 
dye  laser  was  operated  single  longitudinal  mode  at  wavelengths  near  226  nm.  Using  this 
laser,  fully  resolved  absorption  measurements  of  the  NO  lineshape  broadening  parameter  for 
collisions  with  N2  were  obtained.  This  important  data  is  required  in  the  reduction  of  the 
Doppler-shifted  fluorescence  signal.  Extension  of  this  database  to  O2,  H2O,  and  NO 
collisions  as  a  function  of  temperature  will  constitute  an  important  part  of  the  Phase  II  effort. 

Also,  progress  was  made  on  a  high-pressure,  high-temperature  supersonic  burner 
facility.  Nascent  NO  in  the  high-temperature  gases  will  be  supplemented  by  additional  NO, 
mixed  into  the  burner  plenum,  to  generate  a  nominally  Mach  3,  atmospheric  pressure  jet 
flow  with  a  maximum  temperature  on  the  order  of  1500  K.  This  unique  facility  duplicates 
the  fluid  dynamic  environment  of  typical  supersonic/hypersonic  propulsion  test  facilities  and 
will  constitute  a  principle  test  bed  for  the  validation  studies  of  the  Phase  II  portion  of  the 
program.  Hydrogen/air  combustion  was  demonstrated  with  a  maximum  burner  plenum 
pressure  of  80  psi.  The  burner  design  pressure  is  100  psi.  Shakedown  testing  was  plagued 
with  persistent  combustion-driven  pressure  instabilities  which  prohibited  extended  hot  flow 
testing.  Initial  PLIF  images  of  the  NO  distribution  in  a  ~Mach  1.5  cold  air  jet  near 
atmospheric  pressure  were  obtained. 
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ADDENDUM 

Doppler-Shifted  NO  Fluorescence  Measurement 


The  following  section  describes  a  proof-of-concept  experiment  completed  under  the 
Phase  I  program  after  the  submission  deadline  of  the  Phase  I  Final  Report.  The  results  were 
included  in  the  Phase  11  proposal  and  are  described  in  more  detail  here. 


A.i  intn?<iwtiQn 

The  body  of  the  Phase  I  final  report  analyzes  in  detail  the  importance  of  an  accurate 
description  of  the  NO  absorption/fluorescence  lineshape  in  Doppler-shifted  fluorescence 
velocimetry.  In  particular,  the  various  roles  of  pressure  and  temperature  in  determining  the 
width  and  line-center  of  a  static  (zero-velocity)  profile  were  discussed.  The  experiment 
described  in  this  addendum  was  undertaken  to  demonstrate  that,  in  a  strongly  compressible 
flowfield,  a  Doppler-shifted  NO  fluorescence  lineshape  could  be  measured.  Further,  the 
measured  lineshape  could  be  used  to  infer  the  local  gas  velocity  in  the  presence  of  competing 
collisional  and  thermal  broadening  and  shifting  mechanisms.  While  the  Phase  11  program 
will  focus  on  instantaneous  imaging  measurements,  calculation  of  the  gas  velocity  from  the 
NO  PLIF  images  will  require  an  analytical  description  of  the  NO  lineshape  which  is 
demonstrated  in  this  result.  A  technique  similar  to  that  demonstrate  here  will  be  used  in  the 
Phase  n  program  to  obtain  the  baseline  velocity  profiles  for  comparision  with  the  velocity 
imaging  results. 


A.2  Experimental  Description  and  Results 

The  goal  of  this  experiment  was  to  fully  resolve  an  isolated  NO  fluorescence 
lineshape  in  a  compressible  flow  using  the  high  spectral  resolution  ring-dye  laser  described  in 
Section  3  of  the  Final  Report.  The  lineshape  from  the  flowfield  was  compared  to  a  lineshape 
recorded  simultaneously  in  a  static  ceil.  The  frequency  shift  between  the  two  lineshapes,  due 
to  both  the  velocity-dependent  Doppler  effect  and  collisional  shifting,  was  analyzed  to  yield 
the  mean  flow  velocity  at  a  point  in  the  flowfield. 

A  schematic  diagram  of  the  experimental  setup  is  shown  in  Figure  A.I.  The  output 
from  the  frequency-doubled  ring-dye  laser  was  tuned  near  the  02(7)  transition  at  226.72  nm 
and  split  into  two  beam  paths.  One  path  was  directed  through  a  static  cell  filled  with  a 
mixture  of  a  fraction  of  a  percent  of  NO  and  a  few  torr  of  N2.  This  beam  path  was  used  to 
recoi  J  the  fuUy  resolved  absorption  lineshape  as  the  laser  wavelength  was  scanned.  The 
second  beam  path  was  focussed  at  the  centerline  of  a  supersonic  jet  approximately  1  diameter 
downstream  of  the  jet  orifice.  The  laser  beam  was  propagated  at  60  degrees  with  respect  to 
the  flow  axis  to  provide  a  non-zero  scalar  product  of  the  gas  velocity  and  laser  beam  wave 
vector.  Fluorescence  from  the  focal  volume  (~  1  mm^)  was  collected  using  a  spectrally 
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Figure  A.l  Experimental  Setup  for  Doppler-Shifted  Velcscity  Measurement 


filtered  photomultiplier  mbe.  As  the  laser  wavelength  was  scanned,  this  beam  path  generated 
a  Doppler-shifted  lineshape  for  comparision  to  the  static  cell  iineshape. 

The  jet  resevoir  was  an  atmospheric  pressure,  300  K  plenum  of  N2  which  was  seeded 
with  0.5  percent  NO.  The  gas  mixture  expanded  through  the  converging  nozzle  and  emerged 
into  the  10  torr  background  chamber  as  an  underexpanded  Jet.  The  facility  was  exhausted 
through  a  vacuum  pump  in  order  to  maintain  the  10  torr  oackgruund.  At  the  measurement 
location,  the  gas  temperature  is  approximately  70  K  and  the  Mach  number  approximately  4 
(cf..  Reft  11  from  the  Phase  I  Final  Report). 

An  example  data  set  is  shown  in  Figure  A.2,  The  dashed  line  is  the  absorption 
profile  recorded  across  the  static  cell.  The  solid  line  is  the  fluorescence  profile  recorded 
simultaneously  from  the  jet.  The  fluorescence  profile  is  shifted  with  respect  to  the  absorp¬ 
tion  profile  by  -  2.04  G&.  Using  the  flow  temperature  estimate,  the  fluorescence  lineshape 
was  fitted  to  a  Voigt  profile  to  extract  the  local  flow  pressure.  The  best  Voigt  fit  corre¬ 
sponded  to  a  flow  pressure  of  35  torr.  From  Eqs.  (12)  and  (13)  of  the  Final  Report,  the 
collisional  shift  at  these  conditions  was  calculated  to  be  -0.6  GHz.  The  remaining  -1.44  GHz 
shift  is  due  to  the  component  of  the  gas  velocity  vector  aligned  with  the  laser  beam  propaga¬ 
tion  axis.  From  Eq,  (2),  this  velocity  component  is  33,(X)0  cm/s,  corresponding  to  a  directed 
axial  component  of  66,000  cm/s.  At  ~  70  K,  this  corresponds  to  a  flow  Mach  number 
'  3.9,  in  good  agreement  witit  the  expected  value. 

This  experiment  demonstrates  .several  key  features  important  to  the  proposed  Phase  II 
program.  Firstly,  and  most  importantly,  it  demon.strates  that  the  models  developed  in  the 
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Figure  A.2  Comparison  of  NO  Q2(7)  Absorption  and  Fluorescence  Lineshapes  Acquired 
in  a  Static  Cell  and  Supersonic  Jet,  Respectively. 


Phase  I  program  to  describe  the  NO  lineshape  are  accurate.  Extracting  velocity  from 
Doppler-shifted  profiles  in  this  manner  is  critically  dependent  upon  understanding  of 
Doppler-  and  pressure-shifting  as  well  as  Voigt-broadening.  The  measured  shift  was  slightly 
more  than  the  FWHM  of  the  fluorescence  lineshape,  but  nearly  30  percent  was  due  to  the 
pressure  difference  between  the  flow  and  reference  cell.  Without  accurate  knowledge  of  the 
Nj  broadening  and  shifting  parameters,  it  would  not  have  been  possible  to  determine  the 
correct  flow  velocity. 

Secondly,  this  experiment  demonstrated  the  principle  of  the  measurements  which  will 
be  used  to  characterize  the  baseline  velocity,  temperature,  and  pressure  profiles  in  the 
supersonic  burner  experiments  of  the  proposed  Phase  II  program.  Using  a  pulsed-dye  laser 
rather  than  the  ring-dye  laser,  NO  fluorescence  lineshapes  in  the  jet  will  be  recorded  simulta¬ 
neously  with  static  cell  absorption  lineshapes  to  determine  the  local  gas  velocity.  Mean 
temperature  measurements  will  be  made  by  scaiming  the  laser  beam  across  several  isolated 
fluorescence  features  and  comparing  the  measured  intensity  distribution  to  the  Boltzmann 
population  levels  at  a  given  temperature.  The  measured  temperature  profile  will  be  used  to 
fit  the  lineshape  to  a  Voigt  profile  so  as  to  extract  the  local  gas  pressure.  These  baseline 
profiles  will  be  used  to  determine  the  instantaneous  velocity  imaging  measurement  accuracy 
over  a  range  of  velocity,  temperature,  and  pressure  conditions. 

In  summary,  the  results  described  above  verify  the  feasibility  of  our  Phase  II 
approach  by  demonstrating  that  our  understanding  of  the  NO  lineshape  is  sufficiently  detailed 
that  velocity-shifts  may  be  measured  in  the  presence  of  pressure  and  thermal  broadening. 
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